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Abstract 



A detailed description of the tunneling processes within Aharonov-Bohm 
(AB) rings containing two-dimensional quantum dots is presented. We show 
that the electronic propagation through the interferometer is controlled by 
the spectral properties of the embedded dots and by their coupling with the 
ring. The transmittance of the interferometer is computed by the Landauer- 
Biittiker formula. Numerical results are presented for an AB interferometer 
containing two coupled dots. The charging diagrams for a double-dot inter- 
ferometer and the Aharonov Bohm oscillations are obtained, in agreement 
with the recent experimental results of Holleitner et al. [Phys. Rev. Lett. 
87, 256802 (2001)] We identify conditions in which the system shows Fano 
line shapes. The direction of the asymetric tail depends on the capacitive 
coupling and on the magnetic field. We discuss our results in connection with 
the experiments of Kobayashi et al [Phys. Rev. Lett. 88, 256806 (2002)] in 
the case of a single dot. 
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I. INTRODUCTION 



The electronic transport through Aharonov-Bohm rings with embedded quantum dots 
(QD's) is a new subject in mesoscopic physics whose complexity competes with the already 
'classical' problem of persistent currents in closed loops. 

Inserting one dot in a ring Yacoby et al. 1 studied for the first time the transport properties 
of such systems. The observed Aharonov-Bohm (AB) oscillations of the source-drain signal 
as the magnetic field is varied proved that the tunneling current through the dot is partially 
coherent. The experiment presented a striking behavior of the transmittance phase as a 
function of the gate voltage applied on the dot: at each transmittance peak the phase jumps 
by it. Due to the two-lead geometry used in this experiment the conductance obeys the 
Onsager relations and as shown in 2 , this imposes a rigidity of the transmittance phase (0 
or 7r). Later on Shuster et al. 3 ' 4 employed a many lead geometry, the phase evolution being 
obtained as well as the expected Aharonov-Bohm oscillations. The experimental geometry 
was generalized by Holleitner et al. 5 who measured the current through a double dot AB 
interferometer (one QD in each arm of the ring). The main achievement of their setup is that 
the dots can be coherently coupled and hence the transport becomes more complex. They 
have also found AB oscillations of the current and emphasized the formation of coherent 
molecular states in the two dots. Finally, a recent experiment 6 with two-dots AB ring was 
performed in a four-lead geometry, the measured transmittance showing peaks in several 
regimes of the capacitive coupling of the ring. Notably, the phase of the transmittance 
presents the same increment with n when one dot is set to resonance and the capacitive 
coupling of the second dot is varied around a peak. 

A closely related problem is the Fano feature of AB interferometers. As reported in 
Ref. 7 a one dot AB interferometer shows asymmetric line shapes for the transmittance as 
a function of the plunger gate voltage, the typical proof of the Fano effect 8-10 , namely the 
interference between states belonging to continuous and discrete spectra. 

The transport properties of AB interferometers containing QD's were theoretically stud- 
ied by two techniques, the scattering approach and the Keldysh formalism in the tight 
binding picture. The scattering theory was successfully used in 11-13 to describe the physics 
of one-dot interferometers, including specific properties of the transmittance phase. The 
S matrix of the scattering problem is computed by writing the Born expansion for the 
T— operator, the conductance being thereafter obtained via the Landauer-Buttiker formula. 

In the tunneling picture the net current from one lead to another is computed by per- 
turbation theory and non-equilibrium Green function techniques. Within this approach one 
can discuss in detail the co-tunneling spin-dependent processes and finite bias transport 14 . 
As discussed recently in Ref. 15 both approaches are equivalent, in spite of the differences 
between the Hamiltonians (in the tunneling picture the coupling between the ring and dots 
does not appear explicitly). 

The way in which the experiments with AB interferometers can indeed provide the 
transmittance phase is a subtle point that involve the explicit geometry of the leads used to 
break the unitarity of the two-lead system (see 16-19 ). 

In the present work we study systematically the tunneling and coherence properties of 
AB interferometers with QD's, particular attention being payed to the geometry used in 
Ref. 5 . The idea behind the calculations presented below is the following. The transmittance 
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of the interferometer as a whole is first related to its Green function, by the Landauer- 
Biittiker formula. Secondly, it is shown that this Green function can be expressed in terms 
of two Green functions that describe separately the ring and the dots system in the presence 
of the leads. The lead-ring, lead-dots and ring-dot couplings appear as non-hermitian self- 
energies of an effective Hamiltonian. The latter is obtained by the Feschbach formula 20 ' 21 
which is a useful tool when dealing with Hamiltonians of coupled subsystems. We point 
out that this step is necessary in order to obtain detailed information about the complex 
processes within the interferometer. The resonant transport through the device is discussed 
in connection with the spectral properties of the dots system embedded in the interferometer. 
Our approach shows clearly that the important role in the resonant transport processes is 
played by the dots inserted in the ring, the latter providing in turn the suitable geometry 
for quantum coherence. 

We do not consider in this paper the Coulomb repulsion because interaction effects on 
the transport properties of single and coupled dots were studied extensively in the previous 
papers 23-25 and all the analysis presented there remains valid here. The Coulomb interaction 
can be however easily included in our formalism in the Hartree approximation and the 
charging effects are satisfactorily described by this approach (see 11 for a similar discussion 
of the interaction effects in a one-particle approximation). The main topics we consider in 
this work are the tunneling and coherence properties of AB interferometers. The Kondo-type 
effects which are a subject in itself are not discussed here. 

The formalism is presented in Section II. Numerical results are discussed in Section III 
in connection with the experimental findings, a qualitative agreement being found. Since 
we have considered two-dimensional quantum dots the magnetic field dependence of the 
eigenvalues of the coupled dots system is no longer negligible as in the case of a dot modeled 
by a single site. It will turn out that the drift of the levels in magnetic field affects the inter- 
ferometer transport properties. Moreover, the interferometer regime of the device (namely 
the one that exhibits AB oscillations) is more difficult to reveal. Section IV summarize the 
main results and ends the paper. 



II. FORMALISM 

This section contains the theoretical framework we use to study the electronic transport 
in Aharonov-Bohm interferometers with coupled quantum dots. The Hamiltonians are writ- 
ten in the tight-binding (TB) representation which is particularly useful both for describing 
complex geometries and performing numerical computations. We consider a general inter- 
ferometer that consists of an arbitrary number of two-dimensional (coupled) quantum dots 
embedded in a ID mesoscopic ring having N sites. Some of these sites are shared with the 
dots, which are coupled to each other by tunneling Hamiltonians simulating the tunable 
barriers patterned in experiment. The quantum dots are described as finite two-dimensional 
(2D) plaquettes. 

The electrons reach and leave the interferometer through ideal one-channel semi-infinite 
leads attached on the ring or directly on the dots. The Hamiltonian of the whole system 
has the form 

H = H 1 + H L + (1) 
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with 



H 1 = H D + H R + H R °. (2) 

H 1 is the Hamiltonian of the interferometer, H L and H R describe the leads and the trun- 
cated ring, i.e. what is left from it after removing the dots (the notations can be identified as 
well from Fig.l which represents a double dot interferometer). The magnetic flux through 
the ring appears in H R in the Peierls representation as magnetic phases attached to the 
hopping constants along the truncated ring. Their explicit form is obtained by using for ex- 
ample the Landau gauge. H^f n and are the lead-interferometer and ring-dots tunneling 
Hamiltonians: 

#tu 7 n = H LI + H IL = r L J2(\0 a )(a\ + |a>(0 a |), (3) 

a 

H R ° = H DR + H RD = Tj2(e- ilpm \m)(0m\ + e tlpm \0m) (m\) . (4) 

m 

Here t f ,t are the corresponding hopping parameters and Oa (Om) are the nearest sites to 
the contact points a (m) between lead-interferometer and ring-dots. 

ip m is the Peierls phase associated with the pair of sites |0m), \m). Finally H D is the 
Hamiltonian of the coupled dots which is also written in the Peierls representation. It 
includes the individual Hamiltonian of each dot H Dk : 

H D * = -eV k ]T \i)(i\+t D £ e 2 ^'\i){i'\ (5) 

ieQD k <i,i'> 

and the interdot tunneling term H tun (T- mt ), depending on the coupling constant r- mt which is 
the same for each pair of dots {k, k + 1}. We point out that the dots embedded in different 
arms of the ring can be coupled as well, allowing thus complicated electronic trajectories 
within the system. The constant term from the diagonal part of each H° k mimics the 
plunger gate voltages used in experiments to tune the dots to resonance, < > denotes 
the nearest neighbor summation and to is the hopping integral on dots. 

The conductance matrix g a p of a mesoscopic system at zero temperature coupled to leads 
is given by the Landauer-Buttiker formula (see Ref. 22 for a rigorous derivation) 

2 2 

9ap{E F ) = jT a/3 (E F ) = 4jT 4 L sm 2 k\(a\G cS (E F + iO)\f3)\ 2 , a ^ (3, (6) 

where T a p(Ep) is the transmittance, \a), \f3) are sites located on the ring or dots that are 
coupled to the leads, E F = 2t F cosk is the Fermi energy of the leads and t F is the hopping 
integral on leads. The main quantity in Eq. (6) is the effective resolvent of the system in 
the presence of the leads (see Ref. 23 for more details). In our case G c g(z) = (H cS (z) — z)~ l 
where the effective Hamiltonian is defined as 

H cS (z) := H 1 - rlUz) |a r )(a r | + ]T \a d )(a d \\ (7) 

and acts in the Hilbert space of the interferometer Tij only and embodies the influence of 
the leads at the contact points with the ring which we denote {a r } or the dots {ad} through 
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the non-hermitian terms above (these terms represent the so called leads' self energy - see 
Ref. 30 ). The notation (i(z) = (z^f^Jz 2 — 4i|)/2 (=F shows that z belongs to the upper (lower) 
half-plane) and we choose Kez < 2ti- In the sequel we take for simplicity e = h = ti = 1. 

In the previous papers 23 ' 25 we used simpler effective Hamiltonians. In the particular case 
of a single dot weakly coupled to leads (see 23 ) Eq. (6) gives at once the transmittance peaks 
(as a function of the plunger gate voltage V) which are related to spectral properties of the 
dot. Actually the effective Hamiltonian of the dot has resonances with small imaginary part 
located near the eigenvalues of the isolated dot. Similarly, if H c s(z) describes an array of 
identical dots one can obtain and explain the splitting of the Coulomb peaks as a function 
of the interdot coupling r int in terms of the nearly identical spectra of the dots. Moreover, if 
the interdot Coulomb interaction is neglected, the effective Green function can be expressed 
only in terms of one dot Green function by a recursive formula. 

Here formula (6) is not of much use because even if the transmittance peaks can be 
obtained from it by inverting numerically the finite matrix of the effective resolvent, one 
cannot distinguish between the different paths that an electron can follow. Indeed, due to 
the ring geometry and to the coupling between the dots the transport within the device 
is very complex. Besides that, in the experiments the metalic gates defining the dots are 
patterned in the ring arms while the incident electrons from leads enter the ring freely. This 
means that tl ~ 1 thus a discussion in terms of the resonances of H 1 is useless. These 
drawbacks are only apparent and one can rewrite G e ff in a suitable way to recover the 
missing details. The first step is to decompose the Hilbert space of the interferometer as 
Hi — H-d® Hr where 7i D {l-L R ) is the Hilbert space of dots(ring). We denote then by P, Q 
the projectors on these spaces. P, Q are nothing else but families of on-site projections 
from the coupled dots system and the ring. Next, observe that is a small off-diagonal 
perturbation with respect to H D — rld(z) J2a d \&d)(oid\ and H R — rld(z) J2a r \ a r)(^r\ viewed 
as non-hermitian operators in Hd and Hr. This allows us to use the Feschbach formula 20 ' 21 
which expresses the effective resolvent in the following form (see Eq. (6.1) from Section VLB 
in 21 ) 

G eS (z) = G%{z) + (1 - G*(z)QH eS (z)P) ■ (H°(z) - z)' 1 • (1 - PH cS (z)QG R (z)), (8) 
where we denoted G R s (z) := (QH cS (z)Q — z) -1 and the new effective Hamiltonian reads 

H°(z) := PH cS (z)P - PH eS Q(QH eS (z)Q - z)- l QH cS {z)P (9) 



Noticing that in our case PH e s(z)Q = H DR one obtains by straightforward calculations 
explicit formulae for G^ s (z) and H^ s (z) (we use the notation Gij(z) := (i\G(z)\j)) 

G«(z) := (h r - rlUz) Y: \a r )(a r \ - z^j (10) 
H°(z) := H D - rlUz) £ \a d )(a d \ - r 2 £ e-^~^G R mfim/ (z)\ m )(m'\. (11) 

»d m,m' 

The advantage of using the Feschbach formula is that it provides us with two effective 
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resolvents, each one describing individually the pieces that compose the interferometer. G^ s 
describes the truncated ring in the presence of the leads while G® s (z) := (H^ — z)^ 1 is an 
effective resolvent for the embedded system of dots both in the presence of leads and ring. 
We remark that GQ mQm ,(z) (see Eq. (10)) has a nonvanishing imaginary part even if z lies on 
the real axis, due to the non-hermitian coupling to the leads. This happens because Ci(-^) 
is always complex when \E\ < 2ti. By direct computation we express various elements of 
the conductance matrix using Eq. (8) (this time the E F dependence is omitted as well as the 
subscript 'eff' ) 



= Art sin 2 k G* fl „ + r 2 e^^»>Gf ri0m CC r 



_4_2 -2 



= ■■ \t 



R 

a r ,0, 



|2 



g ar ,/3 d = At l t sin k 



9a d ,p d = 4r L sin k 



^a d ,0 d 



(12) 
(13) 
(14) 



In the above equations the summations over m and n are understood. The set of formulae 
(12)-(14) is the main formal result of the paper and the starting point of a detailed discussion 
of the transport processes through the system in terms of the spectral properties of the 
effective Hamiltonian H^. t^ r a is the transmission amplitude from lead a to lead (3 via 
the truncated ring and t®^p r controls the transport via the arm containing the dot(s). In 
the following we consider some particular geometries already used in experiments. 



A. One-dot AB interferometer 

The simplest AB device is realized when there are no leads attached to the dots system 
which in turn is composed of only one dot (this is the geometry used by Yacoby et al. 1 ). Then 
the term containing the sites {ad} vanishes from Eq. (11) and the transport is completely 
described by Eq. (12) that gives the transmittance of the system. 

Let Ei(V) be the i-th eigenvalue of the isolated dot, ipi the corresponding eigenfunction 
and Pi := its associated projection. Note that the eigenvalues Ei(V) and their 

eigenfunctions depend also parametrically on the magnetic field. We describe below 
the resonant transport through Ei(V). The idea is to isolate the resonant contribution in 
the effective resolvent. With the notation P^ := 1 — P i: the effective Hamiltonian can be 
written as 

H° = PiH°Pi + PiH°P^ + PtH^Pi + PtH°Pt, (15) 

and we can apply again the Feschbach lemma for (H^ — z)~ l having PiH^P^ + h.c as a small 
perturbation of PiH&Pi + P^H^P^. Then with the notations G( := (P^H^P^ - z)' 1 
and Di(z) := PiH^P^G^P^H^Pi the effective resolvent reads 

(H°(z) - z)- 1 = Gf + (1 - GiPtH°Pi){PiH°Pi - A - z)~\\ - PiH^Gf) (16) 

and the resonant term is clearly 
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where the resonance width Tj and the shift Aj are flux-dependent quantities, their expres- 
sions being easily identified. Notice also that {^^D^z^ipi) is of order r 4 thus T^z) is of order 
r 2 . Let now z — > E F + iO and suppose that we fix V such that E^V) = Ei — V ~ E F + Aj (E^ 
being the eigenvalue of the dot in the absence of the capacitive coupling). It is clear that the 
main contribution in (16) comes from (PiH^(z)Pi — Di(z) — z)^ 1 since Gf stays bounded 
and the other terms are of 0(t 2 ). The denominator of (PiH^(z)Pi — Di(z) — z)^ 1 reduces 
to resonance width Tj which compensates the multiplicative factor r 2 from the numerator 
of t? D fl • This behavior induces a peak in t„f> and hence in the total transmittance across 
the ring. With these considerations we conclude that for weak ring-dot coupling, whenever 
V comes close to E F — E { for some the transmittance can be written in the form 



= 2 ^ r 8111 kG ^ m G 0n ^ E v _ A _ E iT +0{r). (18) 



Equation (18) is a Breit-Wigner-type formula and gives the transmittance between the 
leads via the quantum dot, as measured in Refs. 1-3. A similar formula was obtained by 
Hackenbroich and Weidenmiiller for a continuous model 12,13 . They supposed that Ei is flux 
independent, which is true only at low magnetic fields and small dots. This assumption 
permits an analytical discussion of the flux- dependence of t®^p r . As we have said, here we 
shall not neglect the effect of the magnetic field on the dot levels. We also point out that the 
'one resonance' form for the dot transmittance was obtained here starting from a many-level 
description of the dot. The rigorous argument for using from the beginning this simplified 
form is that after subtracting the 'resonant' term from the effective resolvent the remainder 
is nonsingular and small. 



B. AB interferometer with a coherent double dot 

When H D describes two coupled dots embedded in different arms of the ring connected 
to two leads (see Fig. 1) we recover the setup of Holleitner et al. 5 In the absence of the 
lead-dot coupling Eqs. (13) and (14) give no contribution thus we are left only with Eq. (12). 
For the simplicity of writing we shall denote Gf r 0m := Gj. and ip m — ip m i := mm /. Since 
Gfjg = in this case the conductance has the form: 



a 4 4 • 2 t 

9ap = 4r L r sin k 



e Lr aa ,U a , b ,U h , 



+ e ie °>'G R a G°,G R p + e ie »">G R n ,GB h G R 



2 



oca' a'b b/3 



(19) 



We remark that the terms G^, b and G^ b , connect points that belong to different dots. The 
effective Hamiltonian in this case is 



H. 



D ( 



J2. \ " 



'G R {z)\m){m'\ 



(20) 



As in the previous section, we are interested in discussing the resonant transport in terms 
of the spectral properties of coupled dots system. Since the double dot Hamiltonian H D 
depends parametrically on the capacitive couplings Vi, V 2 we denote its eigenvalues and 
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eigenfunctions by Ei(Vi, V2) and ipi(V\, V 2 ). The main point is that for suitable pairs {Vi, V2} 
one can bring E^Vi, V 2 ) close to Ej(Vi, V2) for j — This is due to the spectral properties 
of detuned dots. Let us remind here that the detuning consists in applying an additional 
gate potential to one dot while keeping the other gate voltage fixed. In Fig. 2 we show the 
spectrum of the detuned double dot (10x10 sites on each dot) as a function of the detuning 
potential V\ applied on the first dot, for a fixed value of T- mt . For simplicity the undetuned 
dot is not capacitively coupled thus V2 — 0. Obviously, one half of the spectrum shifts 
linearly in V\. The remaining eigenvalues depend weakly on V±, excepting some points of 
avoided crossings. As long as r int 7^ there are no crossings between eigenvalues (on the 
contrary, as shown in Fig. 2 we rather have avoided crossings). Moreover, by perturbation 
theory, near avoided crossings the distance between eigenvalues is of order r- nit . This behavior 
of eigenvalues as functions of V\ and V2 is due to the fact that, roughly speaking, half of 
the eigenvalues 'belong' to QD l5 the other half to QD 2 . As a consequence, when Vi, V 2 are 
tuned such that both Pj(Vi, V 2 ) and Pj(Vl, V 2 ) are near and moreover close to the Fermi level 
we expect that both dots will transmit. Clearly one can study the tunneling through one 
eigenvalue following the same steps as in the analysis of a single dot case. The interesting 
situation is however the one in which the resonant tunneling involves both eigenvalues. In 
the following we show how this appears formally at the level of G D . To this end let us 
introduce the 2-dimensional projection P^ := Pj + Pj, P k being the projection associated 
to the eigenvalue E k (Vx, V 2 ) with k = i,j for % and j fixed. We shall also use the notation 
P± := 1 - Pij. Then PijHgP± + h.c is a perturbation (of C(r 2 )) to PijHgPij + P l jH° s P l j 
and by the Feschbach lemma for G^ s one has 

<?£ = (P£HZP£ - z)- 1 + (H°(z) - z)' 1 + 0{r% (21) 

with 

H°{z) := PijH^Pj - P tJ H°P£(P^H°P£ - z)' 1 P± H°P i3 =: H D (z) - D^z). (22) 

As in the case of a single dot the first term in (21) is small when z — > Ep + ^0 and the gate 
voltages are chosen such that the resonant condition is fulfilled at least around one of the 
two eigenvalues E k (Vi, V 2 ). Then the last step to be done is to write the Dyson expansion 
of (Htft — zy 1 taking PiH^Pj + h.c := V as a perturbation of the 2x2 diagonal matrix 
PiHgPi + PjHgPj: 

GeE — Gf jfiS {z) + Gf jc& (z)VGf jcS + (23) 
where the unperturbed resolvent Gfj eS has the form 



;d , n = \i>j){i>j\ 



The indices i, j were explicitly written for the unperturbed operator (we did not introduce 
another notation). Thus, we have here two resonances of widths Tj, Tj (their expressions are 
complicated but easy to obtain from (22)). It is clear that as long as the dots are coupled 
Ei(Vi, V 2 ) 7^ Ej(V±, V 2 ) thus the two resonances come close but do not cross. Indeed, if for 
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some Vi, V 2 the first term in (24) behaves like l/Tj the denominator of the second term is 
Tint + (A, — Aj) — iTj thus the resonant condition is not strictly achieved. Let us write 
explicitly the second term from the Dyson expansion (23). Since V is off-diagonal one is left 
with: 

|^)(^-|^ + A,W|^)(^l+h.c 
ijM ) ijM ) (Ei(V 1: V 2 ) - A^z) - iT^z) - z) ■ {Ej{V u V 2 ) - Aj(z) - iT 3 {z) - z) ' 

(25) 

Looking at (20) and (22) one notes that the numerator is quadratic in r as well as the widths 
of the resonances Tj, Tj. Thus the perturbative expansion (23) cannot be used in the case of 
decoupled dots since and Ej can cross, the imaginary parts of the resonances are equal 
Tj = Tj = T and Gfj eS VGfj eS behaves also like 1/T. However here we deal with coupled 
dots and as long as r int > r 2 the Dyson series (23) converges and (19) becomes (omitting 
the indexes i,j and eff in Gf^ eS ) 



g aP = Art sin 2 k r 2 (e*-G« G£(?g + e*-" G^G%G^ 

+ e*-*G* G&G*, + e<'< a ,G£ 6 G* + G(r 2 )) f . (26) 

The last formula allows a discussion of the interferometer properties of the device. The first 
two terms represent the direct tunneling through the upper and lower dot, while the terms 
containing G^ b , and describe paths in which the electron tunnels from one dot to the 
other before being transmitted in the leads. At small interdot coupling the cross products 
(a\ipk) {ipkW) 1 (a'\ipk)(ipk\b), {a\ipj){ipj\b) and {a'\ipi){ipi\b') are expected to be small so that 
we can write (keeping only the first two terms from the right side of Eq. (26)) 



9ap = 4r L sin k 




«»G*Mi>i)mb)GK p | e^'G^(a'\^)mb'}G§A | ^ 



E t (V 1} V 2 ) -A t -tT t -z Ej(Vi, V 2 ) - Aj - tT 



(27) 



where TZ collect all the other paths within the interferometer that give smaller contributions. 
Equation (27) will help us to discuss the numerical results from the next section. 

One may notice that in the above analysis the spectral properties of the truncated ring do 
not appear in an essential way in the problem. This could be anticipated from the beginning 
since it is the double dot system that controls the tunneling events. At the formal level, this 
fact is revealed only by using Feschbach formula. 

We mention that our Eqs. (12)-(14) and (18) are similar to the ones obtained previously 
by Hackenbroich and Weidenmuller 11,12 by a scattering theory approach, in the case of a 
single dot embedded in a ring connected to two leads. Here we gave an alternative calculation 
in terms of the Green functions rather than using the S matrix and we generalized the 
discussion beyond the single-dot case. An advantage of our approach is that we do not use 
the Born series which is formally resummed in the scattering approach. 

Let us finally observe that one could not compute the tunneling current through the 
interferometer via rate equation methods used previously in 31 and 32 for weakly coupled 
quantum dots. These approaches would imply in our problem either the computation of the 
probability distributionP(iV, a) characterizing the interferometer in the iV-particle state a, 
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either a perturbative expansion w.r.t the lead-interferometer tunneling Hamiltonian. Since 
the lead-ring coupling constant is rather big the number of electrons in the interferometer 
is not quantized, thus P(N, a) is not well-defined, and the perturbative argument breaks 
down. 

III. RESULTS AND DISCUSSION 

We start this section with the most interesting geometry, the one used by Holleitner 
et al. 5 Following their analysis we first look for the charging diagrams of a ring with two 
identical dots connected to two leads. The dots have 4x5 sites each, while the ring supports 
100 sites. We recall (see also 26 ) that the charging diagrams are plots of the current through 
a system containing two quantum dots as a function of the gate voltages Vi,V 2 applied on 
each dot. In Fig. 3 we present the rhomboids for our system, obtained as follows: for each 
fixed value of V 2 , we varied V\ in the interval shown in the figures and we selected only 
transmittances (i.e. conductances) T12 that are larger than 0.4, which means that what we 
obtain is roughly a map for the peak positions in the plane (Vi, Vi). The magnetic flux is 
fixed. As the interdot coupling increases the diagram changes, due to the usual behavior of 
the transmittance in coupled dots 33 : a regular peak is split into two subpeaks, separated by 
a distance which increases with r int and saturates at perfect coupling (r int = 1). The tunnel 
split peaks of the interferometer transmittance were observed in Ref. 5 both in vanishing 
and strong magnetic fields (see Figs. 4(a) and 4(b) in the cited reference). Figure 4 shows 
our result for the transmittance of the interferometer at uniform capacitive coupling (i.e 
V\ — Vz), fixed magnetic flux and different interdot tunneling constants (here r r and r are 
also fixed). A striking feature is observed in the case of a ring with decoupled dots (the 
dotted line in Fig. 4): the transport is strongly suppressed. This behavior at r- nit = was 
predicted also in Ref. 15 . It differs from the one encountered in the case of double dots 
connected directly to leads, when two subpeaks merge to a single one as r int — > 0. 

Another important aspect of the charging diagram is the drift of the peaks near double 
resonance points, which actually gives the honeycomb pattern. We discuss this in connection 
with Fig. 3(b) using the spectral properties of the detuned dots emphasized in Section II B. 
The traces from the range Vi G (0.11, 0.35) depend weakly on Vi because the corresponding 
eigenvalues of the embedded double dot have this behavior there. A similar behavior is 
observed with the traces in the interval V 1 G (—0.35, —0.11) where the eigenvalues depend 
weakly on V\. This behavior changes drastically when two traces are approaching (around 
point D marked in the figure): they clearly avoid each other, because the eigenvalues of the 
double dot do not cross (r int = 0.2). The avoided crossing is more difficult to discern at 
small interdot coupling, as in Fig. 3(a). The problem of crossing resonances in double-dot 
AB interferometers is discussed in a recent work 29 were it was proved that actually at real 
energies such crossings do not exist. This result coincides with ours. 

In Ref. 5 , the interferometer properties of the system were revealed by the following 
procedure: for a fixed avoided crossing of the charging diagram the current through the 
interferometer was represented as a function of magnetic field. We follow the same strategy, 
by carefully analyzing first what happens to the transmittance at such avoided crossing 
points of the charging diagram. As we have mentioned, the two traces above regions D and 
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C from Fig. 3(b) correspond to two eigenvalues Ei(V\, V 2 ) (i — 1,2) that depend weakly on 
V 2 . Similarly, the traces that approach A and D are associated with Ej{V\, V 2 ). Looking at 
Eqs. (24) and (26) one can notice that as long as V 2 does not align Ej(Vi, V 2 ) to the Fermi 
level, the only terms that produce peaks in the transmittance are the ones involving 
and this happens each time when Ei(Vi, V 2 ) ~ Ep. The main point is that by varying V 2 we 
achieve the resonant condition for the term involving hence both dots will transmit. 

In Fig. 5 we show a detail from the charging diagram in Fig. 3(b), taken in the neighbor- 
hood of almost crossing points A and B. In contrast to the usual picture with sharp peaks 
here we observe (Fig. 5(a)) an asymmetric large tail of the peaks, which shows that in this 
regime the interferometer acts as a Fano system. This happens because one dot (QD 2 ) is 
always set to a resonance thus the corresponding arm of the ring is 'free', providing the 
continuum component for the interference. Formally this is easily understood by looking at 
Eq. (27), because the second term is always large enough and interfere with a quantity (the 
first term) that increases as E^V) approaches the Fermi level. The Fano regime disappears 
quickly as we tune QD 2 away from resonance, the picture of separate peaks being recovered 
(Fig. 5(b)). 

In Fig. 6 a-d the solid lines are plots of the transmittance as a function of V\ when V 2 is set 
close to a resonant value. Remark the sudden drop of the peak after the resonant point and 
the Fano dips. The latter are actually located in the avoided crossing region, which explains 
the small transmittance there. Moreover, the asymmetric tail changes its orientation as V 2 
is slightly varied, i.e the Fano parameter sign changes. Following Kobayashi et al. 7 we shall 
call this feature the electrostatic control of the Fano asymmetric line. In order to explain 
this observation we have to look at the two paths that are involved in the interference. The 
first contribution comes from the resonant tunneling through the upper dot and is given 
in dashed lines in Fig. 6 a-d (the plots were obtained by decoupling the the lower arm of 
the ring from the leads). In this case there is no interference and one gets usual resonant 
peaks. The second contribution is due to the 'background' transmittance of the lower arm 
when V 2 is set close to a resonance and upper arm does not transmit. We illustrate this 
component of transport in Fig. 6e which shows a single peak that appears by varying V 2 
when V\ is far away from resonant values. The points A,B,C,D mark the magnitude of 
the background for four values of V 2 . Clearly, as V\ approaches the resonant points the 
interference becomes possible and the Fano lines appear. By inspecting each of Figs. 6a-d 
in connection with Fig. 6e one gets a description of the line shape for different pairs of 
Vi,V 2 . As long as the transmittance values of the two contributions are located on the same 
side of their corresponding peaks the interference is constructive and the Fano line increase 
up to a maximum which coincides with the resonant peak of the upper arm. In contrast, 
when V\, V 2 are chosen such that the transmittance values are located on different sides of 
the peaks the two path interfere destructively and the Fano line drops to a dip. In particular, 
for V 2 fixed the dips will be located on the same side of the peaks, thus the Fano parameter 
conserve its sign. The appearance of Fano effect in interferometers with embedded dots 
was also discussed in a simple (exactly solvable) model in Ref. 28 , without considering the 
interdot coupling or emphasizing the electrostatic control of the Fano lineshape. 

In the above discussion the magnetic flux was fixed and we have varied V 2 , emphasizing 
the sensitivity of the Fano interference on this parameter. Fig. 7 shows that the shape of 
the Fano line can be equally controlled by varying the magnetic flux, while keeping V 2 fixed. 
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Indeed, as increases from 3.00 to 4.50 the asymmetric tail changes its orientation. This 
effect originates in the field dependence of the dot levels which leads in turn to a shift of the 
background peak. Indeed from Fig. 8a one notices at once that the background peak moves 
to the left as the magnetic flux is varied. In order to make the connection with Fig. 7 we 
marked with points the transmittance values corresponding to the gate voltage V2 = 0.11. 
As a consequence of the magnetic shift the point located at = 3.00 on the left side of 
the peak passed on the upper right side at = 3.80 from where it goes down for = 4.50. 
The same argument used in the discussion of Fig. 6 explains now the change of the Fano 
parameter shown in Fig. 7. 

Fig. 8b shows the - dependence of the resonant eigenvalue of the isolated double dot 
(the line obtained for a vanishing ring-dot coupling, i.e r = 0) and of the eigenvalue of the 
whole interferometer (drawn at r = 0.3). The horizontal lines mark the flux values chosen 
in Fig. 8a. As expected a non vanishing r leads to a hybridization between the spectra of 
the truncated ring a(H R ) and the coupled dots a(H D ). The double dot eigenvalue acquires 
a quasiperiodic modulation with due to the ring geometry. 

By comparing Figs. 8a and 8b we observe that the background peak follows the field 
dependence of the eigenvalue of the isolated double dot and not the one of the interferometer 
eigenvalue. The physical meaning of this behavior is that the resonant transport is controlled 
by the spectral properties of the embedded dots. If the interferometer eigenvalue would 
control the peak position this one should move to the right from = 3.00 to = 3.80, 
according to the trajectory given for r = 0.3. Clearly this is not the case and, up to a shift 
caused by the real part of the resonance the peak obeys the drift of the isolated eigenvalue. 
We stress that this non-trivial effect described above cannot be captured by a theoretical 
model that neglects the spectral properties of the dot in magnetic field. The direction change 
of the asymmetric Fano tail at the variation of magnetic field was experimentally reported 
by Kobayashi et al. 7 in the the case of a one-dot interferometer. We believe that the effect 
we just discussed for the two-dots interferometer is similar. 

We further investigate the behavior of the Fano peaks as a function of the interdot 
coupling. Figure 9(a) shows that the lineshape is very sensitive to this parameter. More 
interesting is the behavior of the interferometer phase along a Fano resonance plotted in 
Fig. 9(b). For weak coupling (and hence for sharp peaks) the phase shows a rapid increase 
by 2-7T. This feature has some connection with the experimental results obtained in a single 
dot interferometer by Kobayashi et al. 7 They reported an increase of 2n for the phase of the 
AB oscillations (we present instead the phase of the transmittance). In our case the second 
dot is set to a resonance so it acts as a free arm of the ring, from where the similarity with 
the one-dot interferometer. By increasing r int the phase becomes a smooth function of V\. 

We now address the problem of AB oscillations. It is clear that they are to be observed 
if both dots are close to resonance, meaning that the gate voltages V\, V2 are suitably tuned 
near some eigenvalues of the double dot. The delicate point is that the eigenvalues depend 
on the magnetic flux through the ring so that for different fluxes one needs different resonant 
values for Vi,V 2 . Otherwise stated, the rhomboids move with ip (not shown). We found that 
for small magnetic fields the changes are not too drastic and that the AB oscillations can 
be captured by monitoring the Fano dip and plotting the transmittance magnitude there as 
a function of the magnetic flux. More precisely, for a given magnetic flux we keep V2 fixed 
and vary V\ in a range that contains only one Fano dip whose transmittance is determined 
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(this is simply the lowest value in the chosen range). Then we repeat the procedure for other 
fluxes, the results being given in Fig. 10. One can recognize at once the Aharonov-Bohm 
oscillations. Their position is slightly shifted due to the phase accumulation within the dots 
(i.e. we express the transmittance as a function of the magnetic flux through the ring while 
the flux encircled by the real trajectories is a bit larger). Notice that the oscillations are 
in phase at all Fano dips. Figure 11 shows that the oscillation amplitude increases as the 
interdot coupling increases. 

We have also investigated a single-dot interferometer (the ring has the same dimension 
while the dot is a 8 x 9 plaquette). When the free arm is decoupled (by making some hopping 
terms zero ) we have the usual peaks corresponding to resonant tunneling via the dot levels 
(Fig. 12(a)). In order to see the Fano features reported by Kobayashi et al. 7 we restore the 
coupling to the arm and we choose the Fermi level such that the 'background' conductance 
of the arm is around 0.3 (if the Fermi level coincides with some eigenvalue of the free arm its 
conductance approaches unity, obscuring thus the contribution of the dot). As expected, the 
symmetric peaks are turned to Fano resonances shown in Fig. 12(b), their correspondence 
being obvious. One notices that the Fano peaks are either wide or very narrow. We have 
checked that this feature remains also valid for other values of the flux and different number 
of sites composing the dot. Remarkably, the Fano parameter takes the same sign between 
succesive peaks. It was suggested recently by Nakanishi et al in Ref. 34 that this feature 
relates to the correlations between the narow and wide peaks. 

IV. CONCLUSIONS 

The main aim of this paper was to present in a unified formalism the basic properties of 
Aharonov-Bohm interferometers with coupled quantum dots. By combining the Landauer- 
Biittiker approach and the Feschbach formula we studied the transport properties of the 
interferometer in terms of the spectral properties of the embedded dots. Our method involves 
only Green functions and can be viewed as an alternative to the scattering theoretical 
approach. In the case of an interferometer with two coupled QD (one QD in each arm of the 
ring) we give a formula (Eq. (27)) which emphasizes the resonant tunneling process through 
a given discrete level from the dots (we recall that along the paper we have considered 
many-level dots). 

Numerical simulations reproduce the stability charging diagrams of two-dot AB inter- 
ferometer reported in the experiments of Holleitner et al. 5 . A careful analysis of the almost 
crossing points of the diagram lead us to several interesting results which are summarized 
in what follows. When the magnetic field is fixed and one dot is set to resonance the inter- 
ferometer transmittance shows Fano lineshapes as a function of the gate voltage applied to 
the other dot. This corroborates with the results of shows clearly the coherent feature of 
the transport through the system. We emphasized and explained the sensitivity of the Fano 
tail to the gate potential on the second dot. 

As we have said, our model includes the effect of the magnetic field on the dot levels. 
It turned out that this effect explains the change of the asymmetric tail as the magnetic 
flux is varied. It would be of great interest to probe experimentally this latter aspect. 
The transmittance assigned to the Fano dips shows Aharonov-Bohm oscillations, in full 



13 



agreement with the observation of Ref. 5 . The influence of the various coupling constants 
was identified. Finally we reproduced the results of Kobayashi et aV . 

The analysis of the 4-lead geometry in view of the very recent results reported in Ref. 6 
is much more complex and requires further investigation. 
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FIGURES 



FIG. 1. Schematic picture of a two-dots Aharonov-Bohm interferometer. The thick solid line 
represents the truncated ring (R). The dashed countour surrounds the interferometer (I), a, (5 are 
the sites where the leads are connected to the interferometer and a, a', b, b' are the contact points 
between ring and dots. 

FIG. 2. Avoided crossings in the spectrum of a 20 x 10 double dot as a function of the detuning 
potential V\ applied on QDi (r = 0.4, $ = 0.15, Tj nt = 0.1). Here is the magnetic flux through 
one cell, in flux quanta. 

FIG. 3. Charging diagrams for the double dot Aharonov-Bohm interferometer (tx = 1, (f> = 3, 
r = 0.3). if is the magnetic flux through the ring, in flux quanta. The traces represent trans- 
mittances bigger than 0.4. (a) r; nt = 0.1, (b) r; nt = 0.2, (c) r; nt = 0.5. The Fermi level is set to 
0. 

FIG. 4. The effects of the interdot coupling r; nt on the electronic transmittance of a double 
dot Aharonov-Bohm interferometer at fixed magnetic flux (p = 3. The same gate potential V is 
applied on each dot tl = 1, r = 0.5. Full line - Tint = 1, long dashed line - Tint = 0.5, dashed line - 
Tint = 0.2, dotted line - T; nt = (transport is strongly suppressed). 

FIG. 5. The structure of the transmittance peaks from Fig. 3(b) around the points of double 
resonance (a). Away from this point one has distinct sharp peaks that turn into Fano peaks at the 
avoided crossing points (b). 

FIG. 6. a)-d) Solid lines: Fano line shapes as a function of V\ for several values of V 2 . The 
Fano tail changes its orientation by passing through a symmetric maxima. Dashed lines: The 
resonant transport through the upper arm of the ring when the lower arm is decoupled from leads. 
Remark the correspondence between the usual peaks and the Fano maxima, (a) V 2 = —0.1175, (b) 
V 2 = -0.1150, (c) V 2 = -0.1125, (d) V 2 = -0.1100 . e) A resonant peak as a function of V 2 . The 
gate potential on QDI was set to V\ = —0.2. The points A,B,C,D correspond to the values of V 2 
chosen in Figs. a)-d). All plots are made for t = 0.3, tl = 1, Tj nt = 0.2, <\> = 3. 

FIG. 7. (Color online). Magnetic control of the Fano interference. As the magnetic field is 
varied the Fano parameter changes sign, t = 0.3, tl = 1, Tj nt = 0.2. 

FIG. 8. (Color online), a) The background peak moves with the magnetic flux; the gate 
potential on QDI was set to V\ = —0.2 . b) The eigenvalue of the decoupled double dot (t = 0.0) 
has a positive slope w.r.t the magnetic flux. The interferometer eigenvalue (t = 0.3) is additionally 
modulated by the hybridization between the truncated ring and the double dot. 
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FIG. 9. The sharpness of the Fano resonances (a) and the phase of the transmittance through 
the interferometer (b), as a function of the interdot coupling: full line - Tint = 0.05, dashed line - 
Tint = 0.15, dotted line - T m t = 0.3. At weak coupling the phase increases rapidly by 2tt while for 
stronger coupling it increases smoothly by 2ir along a Fano resonance. The parameters used are 
V 2 = -0.110, <j> = 3, r = 0.3, t l = 1. 

FIG. 10. The in-phase Aharonov-Bohm oscillations of the transmittance assigned to the Fano 
dips from the region A, B, C, and D in Fig. 3(b). 

FIG. 11. Aharonov-Bohm oscillations in the region D of the charging diagram at different 
interdot couplings: full line - T Ui t = 0.25, dashed line - T mt = 0.2, dotted line - T Ui t = 0.15. Other 
parameters are ti, = 1, r = 0.3, Ep = 0. 

FIG. 12. Transmittance through a single dot interferometer {jl = 1, r = 0.35, = 5, 
Ep = —0.5 ). The dot has 9x8 sites and the ring contains 100 sites, (a) Usual peaks aris- 
ing from the resonant tunneling via the discrete levels of the dot (the free arm of the ring is 
decoupled), (b) The Fano regime: the free arm conducts and interferes with the path along the 
QD. The peaks turn to Fano lineshapes. 
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